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Abstract 



^> ' We study the explicit calculation of the set of superhedging portfolios of contingent claims 

n , . in a discrete-time market model for d assets with proportional transaction costs when the 

underlying probability space is finite. The set of superhedging portfolios can be obtained 
by a recursive construction involving set operations, going backward in the event tree. We 
reformulate the problem as a sequence of linear vector optimization problems and solve it 
O^, by adapting known algorithms. A corresponding superhedging strategy can be obtained 

going forward in the tree. We discuss the selection of a trading strategy from the set of 

f^ , all superhedging trading strategies. Examples are given involving multiple correlated assets 

^ ' and basket options. Furthermore, we relate existing algorithms for the calculation of the 

^D , scalar superhedging price to the set- valued algorithm by a recent duality theory for vector 

^/^ ' optimization problems. 
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^- 
c^ ■ 1 Introduction 

In this paper, a method is provided for an explicit calculation of the set of superhedging portfolios 
of contingent claims in a discrete-time market model for d assets with proportional transaction 
costs when the underlying probability space is finite. 

The set of superhedging portfolios in markets with transaction costs has been characterized 
under appropriate no arbitrage conditions using consistent price systems, the pendant to the 
density process of equivalent martingale measures in markets with transaction costs ([22 ^ 123 1 [35]). 
But until now, an algorithmic approach to calculate the set of superhedging portfolios was only 
provided for d = 2 assets in Roux, Zastawniak [33]. In a recent paper Roux, Zastawniak [3l] 
present an extension of the recursive representation in [33] to d assets, but do not provide 
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an efficient method to implement it and thus lack examples beyond one-period markets. The 
method presented in this paper will overcome this problem. 

Most algorithm so far concerned the calculation of the scalar superhedging price in markets 
with two assets and transaction costs. The scalar superhedging price is the smallest superhedging 
price if a specific asset (usually the domestic currency or a bond) is chosen as the numeraire. 
Often, strong assumptions on the size of the transaction costs or the type of contingent claims 
had to be imposed to provide an algorithm to calculate the scalar replication price and to 
ensure that this price coincides with the smallest superhedging price (see e.g. [3l O [291 EH])- 
Roux, Tokarz, Zastawniak [32], Roux [31] were able to drop those assumptions and develop an 
algorithm to calculate the scalar superhedging price in the two asset case that is based on the 
dual representation of the scalar superhedging price deduced in Jouini, Kallal [21]. The scalar 
superhedging price is important to deduce price bounds. But only the set of superhedging 
portfolios provides full information if one wants to carry out a superhedging strategy starting 
from an initial portfolio that might contain other assets besides the numeraire asset. 

This paper concerns three goals. First, an algorithm is presented to calculate the set of all 
superhedging portfolios in the d asset case, as well as an algorithm to calculate the superhedging 
strategy when starting from an initial portfolio vector in the set of superhedging portfolios. The d 
asset case also allows to consider basket options. Secondly, we will show that the superhedging 
problem in markets with transaction costs leads to a sequence of linear vector optimization 
problems. This generalizes the well known fact that in frictionless markets, superhedging leads 
to a sequence of linear optimization problems. Thirdly, we will show how the above mentioned 
scalar algorithm of [32 ^ [31] can be related to the algorithm presented here by a recent duality 
theory for vector optimization problems. 

To be more precise, we follow the numeraire-free approach of Kabanov |22j . Schachermayer 
[35j to develop an algorithm to calculate the set of all superhedging portfolios. We will show in 
theorem 13.11 that the set of superhedging portfolios SHPt{X) at time t = 0, 1, . . . ,T and node 
uj £ Qf of a claim X can be calculated going backwards in the tree by 

Vw G Qt SHPt{X){u) = X{uj) + Kt{uj) 
VtG{T-l,...,l,0}, Vi^ef^i SHPt{X){uj)= f^ SHPt+i{X){uj) + Kt{u), 

oJGsucc (oj) 

where Kt denotes the solvency cone of the market at time t. The corresponding superhedging 
strategy can be obtained by going forward in the tree. The superhedging strategy might not be 
uniquely determined and different approaches will be discussed on how to choose a strategy. 

The calculation of both, the set of superhedging portfolios and the superhedging strategy for 
a given initial portfolio vector involves the calculation of an intersection and the calculation of a 
sum of unbounded polyhedral sets, which could be realized by methods from computational ge- 
ometry, which are essentially based on the vertex enumeration problem. Instead we reformulate 
the problems as linear vector optimization problems and solve them with Benson's algorithm 
(which also involves vertex enumeration). This reformulation leads to interesting insights into 
the set (or vector) optimization nature of the problem as in each iteration step a solution of the 
dual vector optimization problem is used, where we use a geometric duality approach developed 
in [IE]. 

The first one to connect coherent set-valued risk measures with linear vector optimization 
problems was Hamel p^ , where in section 8 the set- valued average value at risk was formulated 



as a linear vector optimization problem. We pick up on that idea, as the set of superhedging 
portfolios can be seen as a set- valued coherent risk measure (see [15]), but adopt to the fact 
that in our dynamic setting and since superhedging strategies are path dependent, one rather 
formulates the problem as a sequence of linear vector optimization problems in the spirit of 
dynamic programming (see also section [7]) . We obtain an algorithm that grows only polynomial 
as the trading frequency increases when path-independent payoffs are considered. 

It is notable, that even for the determination of the scalar superhedging price the calculation 
of the set of superhedging portfolios at intermediate nodes is necessary. Thus, the calculation 
of the set of superhedging portfolios is not more involved than the calculation of the scalar 
superhedging price. It will turn out that the scalar algorithm based on the dual representation 
of the scalar superhedging price by Jouini, Kallal [21j is related via geometric duality to the 
set-valued algorithm. We will show that the set SHPo{X) can be recovered from the scalar 
algorithm if a certain duality map is applied to a function appearing in the penultimate step of 
the scalar algorithm. 

The paper is structured as follows. Section [2] reviews basic definitions and results about mar- 
ket models with proportional transaction costs and superhedging in those markets. In section [3] 
the above mentioned recursive construction of the set of superhedging portfolios is deduced. 
An algorithm to calculate a superhedging strategy is presented in section 13.11 as well as a dis- 
cussion on how to choose an optimal strategy based on different criteria. Section S] presents 
examples ranging from multi-period binomial models to correlated trees for multiple assets and 
includes European call options as well as basket options. Section [5] reformulates the superhedg- 
ing problem in terms of linear vector optimization. Basic definitions and results from linear 
vector optimization are reviewed in section [6l Section [7] presents the dual representation of the 
scalar superhedging price of Jouini, Kallal [21j generalized to the d asset case, and provides an 
algorithm to calculate the scalar superhedging price for finite probability spaces. This general- 
izes the algorithm of Roux [31], Roux, Tokarz, Zastawniak [32] to d dimensions, which allows 
the consideration of basket options. The relationship between the scalar and the set-valued 
algorithm is deduced in section 17.31 via geometric duality for the corresponding linear vector 
optimization problems. It shows that the scalar algorithm is in fact a set-valued algorithm. 

As a conclusion, in markets with transaction costs, the set of superhedging portfolio vectors 
as introduced in the numeraire-free approach of Kabanov [22], Schachermayer [3^ has to be 
considered at all intermediate steps, even if one is only interested in the scalar superhedging 
price for a chosen numeraire at time zero. In the past, mainly scalar algorithms, i.e., algorithms 
to calculate the smallest superhedging price expressed in a given numeraire, have been studied 
in the literature. In [33], a set- valued algorithm has been obtained in the two-asset case for 
American options and was extended to d assets in [34| . But a lack an efficient implementation 
for d > 2 and thus a lack of examples beyond one-period markets in [M] show that our approach, 
which establishes a connection to linear vector optimization and thus allows the use of well 
established algorithm to solve the problem, is an advantage. In this paper, we present an 
algorithm to calculate the set of all superhedging portfolios that works for any fixed number of 
assets in a discrete time model on a finite probability space. 

The method presented will also prove useful for the calculation of other (dynamic) coherent 
risk measures in markets with transaction costs, price bounds like good deal bounds, or when 
calculating solutions of portfolio optimization problems in markets with transaction costs. These 
problems will be discussed in a separate paper. 



2 Preliminaries 

Consider a financial market where d assets can be traded over finite discrete time t = 0,1, . . . ,T. 
As stochastic base we fix a filtered probability space {^,J^, {J't)'t=oi P)- We assume without loss 
of generality that Tq is trivial. 

A portfolio vector at time t is an J^j-measurable random vector V^ : O — )• IR . The values 
Vt {oj) of portfolio vectors are given in physical units, i.e., the i-th component of the vector Vt {uj) 
is the number of units of the i-th asset in the portfolio at time t for i = 1, . . . ,d. Thus, we do not 
fix a reference asset, like a currency or some other numeraire, and treat all assets symmetrically 
as it was initiated by Kabanov [22]. 

A market model is an {J^t)f=o adapted process {Kt)f^Q of solvency cones. A solvency cone is 
a polyhedral convex cone Kt (w) with IR^ C Kf (u) / IR'^ for all uj e fl and ah t G {0, 1, . . . , T} 
and models the proportional frictions between the assets according to the geometric model 
introduced in Kabanov [22], see also [351 123] . To be more precise, let the terms of trade at time 
t be modeled via an /"j-measurable d x d bid-ask-matrix IIj, as in definition 1.1 in ^35j. That is, 
IIj is a matrix- valued map co i— )■ Hi (a;), denoting the bid and ask prices for the exchange between 
the d assets. The entry tt*-' of lit denotes the number of units of asset i for which an agent can 
buy one unit of asset j at time t, i.e., the pair {~i,'7T^^) denotes the bid- and ask-prices of the 
asset j in terms of the asset i. Furthermore, IIj satisfies 

tt'^>0, l<i,j<d, 
tt" = 1, l<i<d, 

The solvency cone Kt{uj) is spanned by the vectors vr^-'e* — e^ , 1 <i,j < d, and the unit vectors 
e^ , 1 < i < d. Some of those vectors might be redundant. Non-redundant generating vectors 
of solvency cones as well as their positive dual cones in different types of market situations are 
given in the forthcoming paper [27j. Throughout the paper we assume 

IR^\{0}Cinti^i(a;), (2.1) 

which is satisfied if the d assets are liquid in the sense that the exchange rates between any two 
assets are positive and finite. The cone Kt is called solvency cone since it includes precisely those 
portfolios at time t which can be exchanged into portfolios with only non-negative components 
(by trading at the prevailing bid-ask prices at time t). Thus, Kt contains all the information 
about exchange rates and proportional transaction costs at time t, t G {0, 1, . . . ,T}. 

By LJ](J"(,IR ) = L|](r2, J^i, P; IR ) we denote the linear space of all 7"i-measurable random 
vectors X : Q ->■ m!^ for t = 0,1, . . . ,T. We denote by L%Ft,Dt) = L^^{9.,Ft,P;Dt) those 
J^i-measurable random vectors that take P-a.s. values in Dt. 

An IR -valued adapted process {Vtji^Q is called a self-financing portfolio process for the 
market given by {Kt)J^Q if 

Vt G {0, . . . , T} : Vt- Vt-i £ -Kt P-a.s. 

with the convention V-i = 0. 

We denote by At C L[J(J^t,]R ) the set of random vectors Vt'- ^ ^ TR , each being the 
value of a self-financing portfolio process at time T, i.e. At is the set of superhedgeable claims 



starting from initial endowment G IR at time zero. As it easily follows from the definition of 
self- financing portfolio processes, At is a convex cone and 

At = -L^aiJ'o,Ko)-L'dJ'i,Ki) - ... - L^^{Ft,Kt). 

Note that xq + At is the set of self- financing portfolio processes starting from the initial portfolio 
vector y_i = xo G IR at time t = 0. 

The fundamental theorem of asset pricing states that a market given by (^t)^=o satisfies the 
robust no arbitrage condition if and only if there exists a strictly consistent price system {Zt)^^Q 
(see |35] . theorem 1.7). For most of the paper we will assume a finite probability space, then 
the fundainental theorem of asset pricing simplifies: a market given by (-f^t)j=o s^itisfies the no 
arbitrage condition if and only if there exists a consistent price system (.^t)t=o (^^^ I^S E5])- 
The definitions are as follows. The market given by {Kt)^^Q is said to satisfy the no arbitrage 
property (NA) if At n L^{J-t,^^) = {0}. The market is said to satisfy the robust no arbitrage 
property (NA'') if there exists a market process {Kt)f^Q satisfying 

Kt C Kt and Kt\ - Kt Q int Kt P-a.s. (2.2) 

for alH G {0, 1, . . . , T} such that 

lTnL0(JT,]Rt) = {0}, 

where At is generated by the self-financing portfolio processes with 

VtG{0,...,r}: Vt-Vt.i£-kt p-a.s. 

By K^ we denote the set- valued mapping uj i— )■ K^ (uj), where K^ (uj) denotes the positive dual 
cone of the cone Kt (w) for w G fi and t G {0, 1, . . . , T}. Thus, 



K (^) 



ivGWi'^-.yue Kt (w) : v^u > o| 



Let rii^j"^ denote the map uj i— )• rii^^"*' (w) where riK^^ (uj) is the relative interior of Kj^{uj) C IR , 
see e.g. [30]. An IR^-valued adapted process Z = {Zt)t=o ^^ called a (strictly) consistent price 
system for the market model {Kt)t=Q if Z is martingale under P and 

VtG {0,l,...,r}: Zt€ K+\{0} {GiiK+) P-a.s. 

Let us denote the set of all consistent price systems by Z. 

An initial portfolio vector xq G IR allows to superhedge a claim X G Lj](J^r,IR ) if there 
exists a self-financing portfolio process Vt G At such that xq+Vt = X P-a.s., i.e. if AT G xq+At. 
Note that the term 'self-financing' also allows the agents to 'throw away' non-negative quantities 
of the assets, thus the above condition xq + Vt = X P-a.s. really means superhedging and not 
necessarily perfect replication. Let us denote by SHPq{X) the set of all those portfolio vectors 
xo G IR at time t = that allow to superhedge the claim X G L'^{Tt,^ ), i-e., 

SHPo{X) = |xo G IR'^ : a: G xo + At\ . (2.3) 

The set of superhedging portfolios of a claim X G L^{Tt,^ ) can be characterized as follows. 



Theorem 2.1 ( [35j theorem 4.1, [15] corollary 5.4). Under the robust no arbitrage condition 
(NA''), for a claim X G L°(J"t,IR'^) one has 

SHPo{X) = (xo£JR'^:yZ £Z with E[{X'^ Zt)'] < oo : E[X'^ Zt] < x^ZqJ (2.4) 

n {eQ [X] + G (w)) . (2.5) 

{(Q,w)gWl: EQ[X-]<oo} 

The function X M- SHPq{—X) defines a closed set- valued coherent market-compatible risk 
measure as introduced in [15j. Equation (j2.3p is the primal representation of this risk measure 
and corresponds to the very definition of elements xq G IR allowing to superhedge the claim 
X G L^(J^TilR' )• Equation ()2.4p refers to the representation of superhedging portfolios with 
help of consistent price systems as in [22l [231 ES]- Equation (|2.5p corresponds to the dual 
representation of the set of superhedging portfolios as a coherent risk measure (see [15] ) , where 
we denote G{w) = {x G IR : < w'^x}. The dual elements {Q,w) are defined by 

W' = UQ,w)€MlaxM''\{0}:yte{0,l,...,T}:E diag («;) ^| J"* e L'aiTt,K+)\ , 

(2.6) 



and are one-to-one with the set of consistent price systems Z as it was shown in |15j . 

In (j2.6p . ^Af^ = M-i^ {Q,J^t) denotes the set of all vector probability measures with com- 
ponents being absolutely continuous with respect to P, i.e. Qi : J^t — ^ [0, 1] is a probability 
measure on (i7, Tt) such that -^ G L^ (J-r, IR) for i = 1, . . . , d. 

In the following sections we will show that if the underlying probability space is finite, the 
set of superhedging portfolios SHPq{X) as defined in (|2.3p can be obtained by a recursive 
construction, going backward in the event tree, while the corresponding superhedging strategy 
can be obtained going forward in the tree. Based on this recursive structure, an algorithm to 
calculate superhedging portfolios and strategies for finite probability spaces will be developed. 

3 Recursive representation of the set of superhedging portfohos 

Let us assume a finite filtered probability space {^,J-, {J-t)J^Q,P) with P{(jj) > for all cj G J7 
and the usual assumptions regarding the filtration, in particular, {0,^2} = Tq C J^i (1 ... C 
J^T = ^ ■ Let i7j denote the set of atoms of J-^ for t G {0, 1, . . . , T}. 

Let us consider a tree that represents the income of information, i.e., the nodes of the tree 
correspond to the atoms cj G ilj of J-t. In the examples in section |4] we will consider recombining 
trees such that the number of nodes does not grow exponentially with the number of time steps 
to ensure that the algorithm is computationally manageable. However, we do not need this 
assumption now. 

A node w G ^t+\ is called a successor node of a; G ilf (t G {0, . . . , T — 1}) if tD C w. The set 
of successor nodes of cj G Oj is denoted by 

succ (w) = {oj G r^t+i : w C w}. 

Let the market model be described by an adapted stochastic process for the solvency cones 
(-K'i)^o- Oil ^ finite probability space, that is if |0| = A^ G IN, the superhedging portfolios of a 
random payoff X G 1?^(Tt-,^ ) can be written in recursive form. In this setting, we can identify 
L°(Jt,IR'^) withlR'^^. 



Theorem 3.1. If the probability space is finite and the no arbitrage condition (NA) holds true, 
the set of superhedging portfolios SHPq{X) C 1R , defined in (I2.3p . of a claim X G L[](/"t',1R ) 
satisfies 

/ SHPo{X) / IR'^ (3.1) 

and can be obtained recursively via 

yuj^Qr: SHPt{X){uj) = X{uj) + Kt{uj) (3.2) 

yt£{T-l,...,l,0},yuj£nt: SHPt{X){uj)= f| SHPt+i{X){oj) + Kt{uj). (3.3) 

(IiGsucc (w) 

Of course the above theorem is also vaHd for more general probability spaces if one replaces 
the (NA) condition with (NA"^). 

Proof. Let A^o denote the recession cone of a convex set A C ]R , i.e. A^o = {x S IR : Vt > 
A + tx CI A}, see e.g. [30|. The intersection of finitely many non-empty polyhedral convex 
sets Ai, . . . , Am satisfying (^i)oo 2 1R+ for i = 1, . . . , ?7i is non-empty (e.g. the component- wise 
maximum over a selection o* E A, (i = 1, . . . , m) belongs to the intersection). Since Kt{oj) ^ IR_|_ 
for all t G {0, . . . , T} and all a; G O^, the sets SHPt{X){uj) defined by (|3.2p . (|3.3|) are non-empty. 
Let xq G SHPq{X) as defined by the recursive construction (j3.2p . (j3.3p . Equation (j3.3p for t = 
is equivalent to the existence of a Vq G xq—Kq such that Vq G n^en SHPi{X){lu). In particular, 
Vo G SHPi{X){uj) for all w G Jli. Using the definition of SHPi{xj{oj), i.e. dSD for t = 1, there 
exists for each a; G ili a Vi(a;) £ Vq — Ki{uj) such that Vi(a;) G SHP2{X){uj) for all a) G succ (w). 
Going forward in the tree in the same manner and using (|3.2p , one can conclude that there exists 
for each Co G f^T with Co G succ (a;) for some u G ^t-i a Vt(w) G Vt-i{oj) — Kt{Co) such that 
Vt{(jj) = X[Cj). That means, there exists a self-financing trading strategy Vy G xq + Aj- with 
Vt(w) = X{CJ), i.e. X G xo + ^t, which means xq is a superhedging portfolio of X by (12. 3p . For 
the reverse direction, take xq G SffPoC-'^) as defined in (j2.3p . For a given path (wq, • . . ,ut-i) 
with ci;^ G Oj and wt G succ(a;t_i) {t = l,...,r — 1), there exists kt G Kt{u)t) such that 
Xq — /cq — ... — A:t'_i G X{u}) + Kx{Cd) for all Co G O^- with w G succ (ujt-i)- Thus, 

xo-ko- ... - kT-2 G fl SHPT{X){Cd) + A;t-i, 

[DGsucc (tjy_i) 

where SHPt{X) is defined as in (13. 2p . Since this holds for any wt-i G r^T-i with ujt-i £ 
succ (ci;t_2), we have 

xo-ko- ... - A:t-2 e n SHPt{X){Co) + irT-i(w) = 5/7Pt-i(^)(^) 

ljGsucc (w) 

via (j3.3p for all w G VLt-i with w G succ {uot~2)- Succeeding like this reveals xq G SHPo{X) as 
defined by the recursive construction (13.21) . (13. 3p . This proves the equivalence of (12. 3p and the 
recursive definition of SHPo{X) in ()3.2p . ()3.3p . By the fundamental theorem of asset pricing, 
no arbitrage implies the existence of a consistent price system Z G Z (see |24t [35] ) . For a finite 
probability space, £'[(X^Zj')~] < oo is always satisfied and from theorem 12.11 SHPn(X) ^ IR 
follows. D 



Clearly, SHPt{X){oj) is the set of superhedging portfolios of X G L^{Tt, IR ) at time t at node 
cj G rjf, t G {0, 1, ...,T} and (j3.ip is satisfied likewise. As a consequence of theorem 13.11 for all 
t G {0, 1, ...,T} and all uj G fij, SHPt{X){uj) is a non-empty polyhedral convex set satisfying 

int(5i7Pi(X)M)^DlR^\{0}. (3.4) 



The recursive structure ()3.2p . (13. 3p provides a geometric intuition for designing an algorithm 
to calculate the superhedging portfolios. The set of all possible superhedging portfolios can 
be calculated going backwards in the tree. The operations involved are the calculation of in- 
tersections of unbounded polyhedral sets and the sum of such a polyhedral set with a convex 
cone. These operations could be realized by methods from computational geometry, which are 
essentially based on the vertex enumeration problem, see e.g. [ll[6j. For several reasons, which 
will be discussed in section [5l we proceed in a different way. We reformulate in section [5] the 
problem as a sequence of linear vector optimization problems and solve them with Benson's 
algorithm, i.e. in each iteration step a linear vector optimization problem is solved. Benson's 
algorithm also involves vertex enumeration but additionally several linear programs. Note that 
for path independent payoffs, the recursive pricing procedure described in theorem 13.11 grows 
only polynomial as the trading frequency increases. 

It is well known that in the frictionless case, the superhedging problem is a sequence of 
linear (scalar) optimization problems (see e.g. chapter 7.1 in [10 1). Thus, the interpretation of 
the superhedging problem in markets with frictions as a sequence of linear vector optimization 
problems makes very clear how the structure of the problem changes if frictions are allowed and 
how the frictionless case fits into this framework. 

Note that the same geometric intuition as in ()3.2p . (j3.3p appears if one is only interested 
in calculating the scalar superhedging price (i.e. the smallest superhedging price in a given 
currency or numeraire) when transaction costs are present. Even in the scalar case, it cannot 
be avoided to use the set-valued operations (intersection and sum of polyhedral sets). As a 
consequence of (|3.4p . the polyhedral convex sets SHPt{X){u}) can be expressed as the epigraphs 
of piecewise linear functions / : IR ~ — )• IR. For example, in the two asset binomial model, this 
recursive structure can be rediscovered in the sequential problem Qt, p. 71 in [3] 

On the other hand, for the purpose of comparing our method with algorithms developed for 
calculating the scalar superhedging price based on a dual description as in [32l |33l [31] it is quite 
helpful to reformulate ()3.2p . (13. 3p as a sequence of linear vector optimization problems. As we 
will see in detail in section [5l in each iteration step a solution of the dual vector optimization 
problem is used, where we use the geometric duality approach developed in jl8]. As it will 
turn out, this geometric duality provides a link between (j3.2p . (j3.3p and above mentioned scalar 
algorithms and allows to recover the whole set of superhedging portfolios from intermediate 
results of the scalar algorithm. It is a somewhat surprising insight that the calculation of the set 
of superhedging portfolios is not more difficult than calculating the scalar superhedging price 
in markets with transaction costs. The scalar problem will be made precise in section [7] and its 
relation to theorem 13. II will be exploit in detail in section [7l3l 

Remark 3.2. Note that it is sufficient to develop an algorithm for superhedging portfolios 
and superhedging strategies since the set of subhedging portfolios SubHPQ{X) of a claim X G 
L^{J^T, IR ) is just the negative of the set of superhedging portfolios of —X 

SubHPo{X) = -SHPo{-X). 

The corresponding strategy for the buyer of a claim X is to superhedge —X. 



Example 3.3. Let us consider a simple introductory example: a one period binomial model 
with non-constant proportional transaction costs, where the set of superhedging portfolios has 
multiple vertices. We will use this example to illustrate the algorithm and to explain differences 
between the scalar and the set- valued approach. Note that the transaction costs are chosen to 
be quite large, just for the purpose of obtaining illustrative pictures. 

Let asset be a riskless cash account and let us assume for simplicity that interest rates are 
zero. Asset 1 is a risky stock, whose bid-ask prices [S^, Sf) at time t = and t = T are modeled 
as follows: 

(20, 26) 
(18, 25) 

(16, 23). 

We consider a digital option, more specifically an asset or nothing call option with physical 
delivery and strike K = 24. The payoff is given by 

X (u) = (Xi {u;) , X2 {u)f = (0, /{5„ >^| {u;)f. 

Thus, the payoff in the up-node is X{uji) = (0,1)"^ and in the down node X{uj2) = (0,0)"^. 
The calculation of the set of superhedging portfolios by the recursive procedure described in 
theorem 13.11 and illustrated in figure [1] reveals that SHPq{X) has two vertices, one at (0,1)-^ 
and one at (—80, 5) and a recession cone equal to the solvency cone Kq at initial time which 
is generated by (—18,1) and (25,-1) . Figure [2] shows the set of sub- and superhedging 
portfolios. 

The scalar superhedging price is given by 25$ and corresponds to the buy and hold strategy 
that superreplicates the claim. The strategy is to transfer the initial position (25, 0) into the 
vertex (0, 1)"^ of SHPq{X) at initial time and hold this portfolio until terminal time. 

However, the knowledge of the scalar superhedging price and the corresponding strategy 
does not give information about optimal strategies if one already owns some shares of the stock. 
For example, if one owns 5 stocks and is short 80 units cash at initial time (the portfolio 
corresponding to the second vertex of SHPq{X)), one cannot reach the scalar superhedging 
price, that is the portfolio (25, 0) , by selling the stock at initial time. This is illustrated in 
figure El Therefore, the portfolio (—80, 5) does not allow to superreplicate if initial positions 
are "cash only" positions. 

Consequently, if initial positions in several eligible assets (here stock and cash) are allowed 
instead of only one (like cash) the cost of superreplication can be reduced: The portfolio 
(—80,5)"^ G SHPq{X) clearly allows to superreplicate (in this example even to replicate) the 
claim. The scalar superhedging price gives information about price bounds, but for the purpose 
of actually carrying out a superhedging strategy, only the set of all superhedging portfolios gives 
full information on optimal strategies. 

3.1 Calculation of the superhedging strategy 

From the proof of theorem 13. H one can see that the superhedging strategy, when starting from a 
particular element in the set of superhedging portfolios, can be calculated going forward in the 
tree. In the presence of transaction costs optimal superhedging strategies are, in general, path- 
dependent even for path-independent payoffs. But one only needs to compute the superhedging 
strategy along the realized path, one step at a time, in real time. 
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Figure 1: Example I3.3t Illustration of the recursive algorithm (j3.2p . (jS.Sp of theorem 13. 1[ 
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Figure 2: Example 13. 3t The set —SHPq[—X) of subhedging portfolios and the set SHPq{X) of 
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Figure 3: Example 13.31 Pi can be exchanged into P4 = (10,0)^, but not into the scalar super- 
hedging price P3 = (25,0)^. 
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Remark 3.4. A self-financing trading strategy starting from an initial portfolio vector xq G IR 
is a sequence ko,ki, ...,kT of trades at time t = 0,1, ... ,T with kt £ Kt for all t. The j-th 
component of kt gives the number of shares of asset j to be bought/sold at time t and paying 
the prevailing transaction costs as given by Kt- Of course, such a self-financing trading strategy 
is one-to-one with a self-financing portfolio process Vo,Vi, ...,Vt via 

Vo = xq- ko and Vi G {1, . . . , T} : Vt = Vt^i - h 

as it can be seen from the definition of a self-financing portfolio process. 

Theorem 3.5. Let the probability space be finite. For a claim X G L[](J^t,1R ) and a given 
initial portfolio vector xq G IR with xq G SHPo{X), there exists a superhedging strategy for a 
path {ujQ, . . . ,ujt) with ut G 0,t and ojt G succ {oJt-i) (i = 1, • • • ,^)- Such a strategy is given by 
a self-financing portfolio process Vq,Vi,...,Vt defined by 

Vb G ({xo} - A-q) n fl SHPi{X){uj), (3.5) 

(DGsucc (oJo) 

Vt£{{Vt-i}-Kt{uJt))r\ fl SHPt+i{X){Cj), (3.6) 



lilSSUCC (cJt) 



for allt G {!,..., T}. 



Proof. The assertions follow from the recursion (j3.2p , (|3.3p and the proof of theorem 13.11 In 
particular, it follows that the sets in (13. Sp . (13. 6p are non-empty and thus the existence of a 
superhedging strategy follows. D 

Note that a superhedging strategy is typically not uniquely determined. In the following we 
will discuss how to choose an optimal superhedging strategy based on specific criteria. Since 
the strategies are in general superhedging and not replication strategies, it might be possible 
in certain scenarios to withdraw cash or assets without endangering the superhedging criteria. 
Thus, one possible criterion to choose a strategy might be to withdraw as much as possible of 
certain assets at intermediate points in time. A different criterion would be to avoid trading if 
possible. It is also possible to explore the tradeoffs between multiple objectives in a multiple- 
criteria strategy selection process using vector optimization. 

At each time t = 0, 1, ..., T, choosing a superhedging strategy at a node cj G 0^ with endow- 
ment X = Vtiyj) G SHPt{X){(jj) means choosing a triplet {v,y, k), where v = Vt+i G SHPt+i is 
the portfolio after any trades k G Kt{Lo) and withdrawals of a portfolio y G IR at time t. That 
means, Vj+i = Vt — y — k and the portfolio V^+i is hold from time t to t -|- 1 and presents the 
initial endowment (before trades are made) in the next iteration step. 

Max-^vithdrawal superhedging strategy. 

Let us assume an investor following a superhedging strategy wants to withdraw as much of a 
certain portfolio y G IR_|_ \ {0} as possible at each intermediate point in time. 

Let i G {1, . . . ,T}, oj £ 0,t, X £ SHPt{X){ijj) and let Kt{Lo) G IR"'^'' be a matrix containing 
the s generating vectors of Kt{uj). We assume that inequality representations of the polyhedral 
sets SHPt^i{X)(uj) are known. An algorithm to compute them will be given in section[5j Solving 
the following LP with variables (f,a,z) G IR x IR x IR^, we obtain a portfolio v according to 
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theorem 13.51 which has the property that a maximal amount of the portfoho y is withdrawn at 
time t. lit < T, we take 

maxa s.t. ve Q SHPt+i{X){Cj), v + ay + Kt{(jj)z = x, a>Q ,z >(). (3.7) 

aJGsucc ibj) 



At time T, we solve 

maxa s.t. v > X{(jj), v + ay + Kt{uj)z = x, a > , z > 0. 



(3.8) 



From a solution {v,a,z), we get the new portfolio Vt+i = v. The portfolio y = ay ^ IR,^ 
describes the withdrawal and k = Kt{uj)z is the corresponding trade. 

Theorem 3.6. If the market satisfies the no arbitrage property (NA), then there exists a solution 
{v,a,z) for dSZl) and (ISlHD . 

Proof. The feasible set is non-empty by theorem 13.51 and the fact that any nonnegative with- 
drawal y belongs to Kt{uj). Assume there does not exist a solution {v,a,z) for (j3.7p or (|3.8p . 
i.e. the value of the problem is -|-oo. Denoting by C = {SHPt{X){u;))oo the recession cone of 
SHPt{X){uj), we obtain y £ -C. But y € int C by dS^D- Hence G int C and thus C = IR'^. It 
follows that SHPt{X){uj) = IR'^, which contradicts theorem EH D 

An important special case is the max-cash superhedging strategy, which is obtained by setting 
y = (1,0, .. . ,0)"^", where the first component is assumed to correspond to the cash account of 
interest. 



Min-trading superhedging strategy. 

Another possible strategy to choose a superhedging portfolio is avoid trading whenever possible. 
The reason to minimize trading activities is that trading in markets with transaction costs is an 
irreversible process. The natural question how to measure trading activities meaningful leads 
to the problem to weight trading in case of multiple assets with different bid-ask spreads. For 
the case the generating vectors of the solvency cone are given by 



K, 



( {stf 



\ 



\s\Y 



{St? 

1 

-1 



<s\f 



1 



{stY-^ 






\s\) 






"'\ 



(3.9) 



where (5^)*, (5" )* are the bid and ask prices -n- and vr of asset 



I, I 



2, ...,d in units of asset 



1, we will assume that the weights are chosen proportional to the corresponding bid-ask spreads 
7^ := {Sfy — [S^y. Solvency cones of the above form appear for example is all assets are denoted 
in the same currency and the riskless asset is the currency itself or a bond with zero bid-ask 
spread (see [27] for more details). We obtain a weight vector 



{7,^7.^7.^7.^...,7^^7^^y 



(3.10) 
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For t G {1, . . . , r} and uj ^ 0,t and a portfolio x E SHPt(X){ij) we solve one of the following 
LPs. If i < T, we take 

min7^2 s.t. ve f] SHPt+i{X){0), v + Kt{uj)z = x, z > 0, (3.11) 

oigsucc (a;) 

where we denote by Kt{oj) the generating vectors of Kt{u)). At time T we solve 

niin7 z s.t. v > X{uj), v + Kt{oj)z = x, z > 0. 

Since the feasible sets are non-empty by theorem 13.51 and the objective function is bounded 
below by zero, there always exist solutions to (I3.10p and ()3.1ip . 

Bi-criteria strategy selection. 

Naturally there is a tradeoff between a max-withdrawal and a min-trading strategy, which can 
be displayed to the agents by solving a vector optimization problem with two objectives. 

Using the above notation we solve in case of t < T the following linear vector optimization 
problem: 

min(-a,7^z)'^ s.t. ve f] SHPt+i{X){0), v + ay + kt{uj)z = x, a > 0, z > 0. {3.12) 

aJGsucc (lu) 

At time T we solve 

min(— a, 7 z) s.t. v > X{uj), v + ay + Kt{ijj)z = x, a > 0, z > 0. (3.13) 

If the market satisfies the no arbitrage property (NA), the problems (I3.12P and (]3.13p are 
feasible and bounded in the sense that there exists I S IR such that {—a, 'y'^z)'^ S {/}+IR_,_ for all 
feasible points {v,a,z) (compare theorem 13. 6p . From |26i theorem 5.20] and the considerations 
in section 5.1 in [2S] we conclude that a (finitely generated) solution always exists. Solutions to 
()3.12p and (j3.13p can be obtained by Benson's algorithm, compare \26\ chapters 4 and 5]. The 
solution concept is also outlined below in section [H 

An example involving the max-cash, min-trade and the bi-criteria strategy selection is given 
in example 14.41 in the next section. 

4 Examples 

If one is interested in the superhedging portfolios for initial positions in only a few of the d 
assets, for example in just a few currencies, or even in just one currency, one would calculate 

SHP^^X) = SHPo{X) n M, (4.1) 

for M = (X^jgj Sis', Si G IR}, where the assets i £ I with I C {1, ..., d} are the ones of interest. 
This simply involves one more operation (the intersection) in the algorithm. The introduction 
of M is quite useful if the number of assets d is very high, one owns just a few of the d assets 
and thus is interested in superhedging portfolios just starting with those assets, or if SHPq^X) 
is too complex to be visualized. In particular, the i-th component of the vertex of SHPq {X) 
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for M = {se*,s € IR} coincides with the smallest super hedging prices if asset i S {l,...,d} is 
chosen as the numeraire, thus coinciding with the scalar superhedging price 7rf{X). 

A second possibility to calculate the scalar superhedging price vr?(X) w.r.t. a numeraire 
asset i is by normalizing the inequality representation SHPo{X) = {x € IR | Bx > b} obtained 
by the SHP- Algorithm in section O below in the following way: Transform Bx > b into Bx > b 
such that each element in the ith row of B is equal to 1. This is always possible by (|3.4p and 
SHPq{X) 7^ IR . Then, the largest component of b is 7rf{X) w.r.t. the numeraire asset i. This 
relation will also play a role when the connection between the set-valued approach and scalar 
algorithms is established and will be discussed in more detail in section 17.31 

Both methods provide a way to compare our results with the algorithm for the d = 2 case 
of j32t [3T| for calculating the scalar superhedging price when the numeraire asset is the riskless 
asset (see example 14. 2p . 

4.1 Two asset case 

Example 4.1. Let us consider a digital option similar to example 13. 3| but in a multi-period 
framework and smaller transaction costs. Let asset 1 be a riskless bond B with an annual 
interest rate of 3%, face value Bt = 1, maturity one year, frequent compounding with n = 100 
time intervals, i.e. Bq = {l-\- -)"" and no transaction costs for the bond, i.e. B^ = Bf = Bt for 
all t. Let the mid-market stock price S follow a Cox-Ross-Rubinstein binomial model, 

St = EtSt^l, 

for t = 1,...,T, where £i,e2,--- is a sequence of independent identically distributed random 
variables taking two possible values e'^ or e~'^ , where At is the length of one time step. The 
initial stock price is 5*0 = 18, volatility a = 0.2, and transaction costs are constant A = 0.04. 
Let the bid and the ask price at time t, respectively, be given by 

St = St{l-X), S^ = St{l + X). (4.2) 

An asset or nothing call option with physical delivery, maturity 1 year, strike K = 19 and payoff 

X (a;) = (Xi {u;),X2 {u;)f = {o,I{sa>K} (^))^ 

is considered. The set SHPq{X) has two vertices, one at (0, 1) and one at (—24.92, 2.39) and a 
recession cone equal to the solvency cone Kq at initial time which is generated by [—Sq/Bq, l) 

and (^oV^o, -l)'^. 

The scalar superhedging price in the numeraire asset (the bond) is 7Tq{X) = 19.29 units 
of the bond, which corresponds to a scalar superhedging price in the domestic currency of 
7r"(X) = 18.72 = ^o and the corresponding strategy is the buy and hold strategy. Note that 
in contrast to the trivial strategy one obtains for the scalar superhedging, the superhedging 
strategy can be more involved, when the initial portfolio vector is not cash-only. 

Example 4.2. Let asset 1 be a riskless bond B with an effective interest rate of r^ = 10%, 
frequent compounding, face value Bt = 1, maturity 1 year, i.e. Bq = {1 + re)~^, and no 
transaction costs for the bond, i.e. B^ = Bf = Bt for all t. Let the stock price S follow a 
Cox-Ross-Rubinstein binomial model as in example 14.11 The initial stock price is Sq = 100, 
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volatility a = 0.2, maturity T = 1 year and transaction costs are constant A = 0.00125. Let the 
bid and ask prices at time t be given as in ()4.2p . Consider a call option with maturity 1 year, 
physical delivery and strike K = 80 whose payoff is a function on the mid-market price, i.e. 

X iu) = (Xi (u) , X2 {u)f = {-KI{Sr>K} (u;) , I{Sr>K} {^)Y • 

The set of superhedging and subhedging portfolios is given by its vertices and recession cones. 
For different values of n, the vertices (in units of bond and stock) are recorded in table [TJ 
The recession cone of SHPo{X) is always Kq, generated by (— 5g/-Bo,l) and (Sq/Bq,—!) , 
whereas the recession cone of —SHPq{—X) is —Kq. The scalar price bounds ir^ {X) , ir'^ (X) in 
the domestic currency are also recorded in table [H 

For comparison purpose, we also give the scalar price bounds n'' (X) , ir'^ {X) if there are no 
transaction costs at i = as considered in [32l [HI [U [28] . We are able to replicate the scalar 
results as given in table 1 of [32] and table 3.1 and 3.2 of [31], where the different values of 
the parameters K and A are K E {80, 90, 100, 110, 120} and A G {0%, 0.125%, 0.5%, 0.75%, 2%}. 
Minor deviations (all less than 0.001) from table 3.2 of [31] appear in a few instances for the bid 
prices in the n = 1000 case and one deviation of 0.014 for n = 1000 that is recorded in table [H 
The case n = 1800 was not considered in [32[ I31j . Here, we just present the results for i^ = 80 
and A = 0.125%. 



A = 0.125% for all t | 


n 


6 


13 


52 


250 


1000 


1800 


vertex of -SHPo{-X) 


/ -74.434 \ 
V 0.953 J 


/ -74.699 \ 
V 0.956 J 


/ -75.477 \ 
V 0.962 J 


/ -76.348 \ 
V 0.969 J 


/ -78.049 \ 
V 0.983 J 


/ -79.049 \ 
V 0.992 J 


lower price bound n''{X) 


27.552 


27.537 


27.462 


27.381 


27.249 


27.191 


vertex of SHPo{X) 


/ -73.814 \ 
V 0.948 J 


/ -73.857 \ 
V 0.949 J 


/ -73.857 \ 
V 0.949 J 


/ -72.856 \ 
V 0.941 J 


/ -71.244 \ 
V 0.929 J 


/ -70.209 \ 
V 0.921 J 


upper price bound 7r°(X) 


27.854 


27.866 


27.872 


27.994 


28.213 


28.370 


A 


= 0.125% for t = 


1,...,T, but no tr 


ansaction costs at 


t = asin [32l[3ll[l[28! 




n 


6 


13 


52 


250 


1000 


1800 


lower price bound tt''{X) 


27.671 


27.656 


27.582 


27.502 


27.372" 


27.315'' 


upper price bound ■7t'^{X) 


27.735 


27.747 


27.753 


27.876 


28.097 


28.255'' 


" differs from value 27.386 


in [3i] 










*> not considered in [32;, :3T] 













Table 1: set- valued and scalar sub- and superhedging portfolios of European call options 

Note that, if the bond is chosen as the numeraire asset, the scalar superhedging price vr"(X) 
is given in units of the bond, and one needs to multiply it by Bq to obtain the scalar superhedging 
price vr"(X) in the domestic currency that is recorded in table [TJ It is worth pointing out that 
there are parameter constellations that lead to multiple vertices for the set of superhedging or 
subhedging portfolios. For example, —SHPq{—X) for A = 2%, K = 110 and n = 52 has 8 
vertices given by the columns of the following matrix 



34.743 


-48.097 


-79.757 


-88.323 


-91.778 


-84.331 


-54.520 


-41.461 


0.322 


0.445 


0.732 


0.809 


0.840 


0.774 


0.504 


0.384 
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with a scalar subhedging price of vr (X) = —0.023 (in the domestic currency) if transaction costs 
are considered at all time points, and 7r''(X) = 0.865 if no transaction costs are considered at 
t = 0. The set -SHPo{-X) for A = 2%, K = 110 and n = 250 has 3 vertices given by 

2.370 -107.125 -110.107 
-0.036 0.974 1.001 

with a scalar subhedging price of vr^(X) = —1.546 if transaction costs are considered at all time 
points, and 7r''(X) = —0.038 if no transaction costs are considered at t = 0. Note that negative 
bid prices might occur when physical delivery is considered in markets with transaction costs. 
This issue was discussed and resolved in [29], see also remark 3.30 in |31j . 

4.2 Multiple correlated assets and basket options 

We are interested in the set of superhedging portfolios of options involving multiple correlated 
assets. We will use a multi-dimensional tree that approximates ad— 1-dimensional Black Scholes 
model for d — 1 risky assets, where the stock price dynamics under the risk neutral measure Q 
are given by 

dSf = Siirdt + aidW^), i = 1, ...,d - 1 

for Brownian motions W^ and W^ with correlation pij £ [— Ijl] for i ^ j- We will follow the 
method in [25] to set up a tree for the correlated risky assets by transforming the stock price 
process S into a process Y with independent components. This tree will have 2'^~^ branches in 
each node and will be recombining with [t + 1) nodes at time t with t G {0, 1, ..., T}. Thus, 
a node can be identified by an index (t, ji, ■■■,jd-i) for t G {0, 1, ...,r} and 1 < ji < t + 1 for 
all i £ {1, ...,d — 1}. For d = 3 the nodes at time t can be described by the indices (ii,j2) of 
the elements of a matrix Mt G IR"^ )^(*+ ). Xhe values of the process Y at such a node can be 
obtained as follows. 

Let S be the covariance matrix of the log asset prices and GG^ = S be the Cholesky 
decomposition of S. Let n be the number of time intervals and At the length of one time interval. 
Let us denote a = G^^{r — ^o"^). The initial value of the process Y is given by Yq = G^^{Xq) 
with Xq = (log(5'Q), ...,log(S'Q~ )). The value of the process Y at node (t,ji, ...,jd~i) is given 

by 

y/ = yj + toiAt + (2ji - t - 2)\/At, i = l,...,(i-l. (4.3) 

for t G {0,l,...,r} and I < ji < t + I. We omit the index (ji, ..., jd-i) for Yt{ji, ...,jd-i) 
and hope not to cause confusion. The value of the original stock price vector 5 at this node 
{t,ji,...,3d-i) is 

5i* = exp(^^*)\ i = l,...,d-l. (4.4) 

Now, let us assume for simplicity that the proportional transaction costs are constant for each 
of the risky assets and are given by A = (A^, ..., A"'^^). Thus the bid and ask prices at node 
{t,ji,-,jd-i) are given by 

{Str = Si{l-\') {S'iy = Sl{l + X'\ i = l,...,d-l. (4.5) 
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Furthermore, let us assume there is a riskless asset with dynamics {Bt)f^Q. Transaction costs in 
the riskless asset can be incorporated by considering bid-ask prices B^ < Bf for i = 0, 1, ...,T. 
For d = 3, if both risky assets are denoted in the domestic currency (the currency of the riskless 
asset), if A^, A^ > 0, and if we assume an exchange between the two risky assets can not be made 
directly, only via cash by selling one asset and buying the other, we obtained the following tree 
model for the solvency cone process Kt- At node {t,ji,J2) the generating vectors of the solvency 
cone are given by vr'-'e* — e^ , < i,j < 2 (see section [2]), i.e., by the columns of the following 
matrix 



is?r {s^y {Sir (s!? 

B\ Bf B'' 



^00 
-110 1-1 

y u u i i j^ ^^,^2 y 



(4.6) 



If there are no transaction costs for the riskless asset, i.e. B^ = Bf, the last two generating 
vectors in (j4.6p are redundant and can be omitted. Non-redundant generating vectors of positive 
dual cones of solvency cones Kt are given in [27] ■ 

Note, if there is an assets denoted in a currency different from the domestic currency, one 
needs to model the exchange rate between domestic and foreign currency as well and obtains 
a model with one more risky asset. The generating vectors of the solvency cones in higher 
dimension and without the above assumptions can be obtained analogously, see definition 1.2 in 

m- 

Also, if one of the risky assets is a currency, one would rather use a discrete approximation 
of a mean-reverting process than a geometric Brownian motion for this asset. In this case 
the model for the process Y with independent components would be similar to above, but the 
transformation ()4.4p needs to be adapted to the new setting, see for example |19j . 

Random proportional transaction costs for an asset can be modeled analogously, by treating 
them as another (correlated) risky asset. In ()4.5p one would just replace the constant A by the 
value of the stochastic process A at node (t, ji, •••, jd-i)- 

Example 4.3. (exchange option) Let us consider a European option in which at expiration, the 
holder can exchange one unit of asset 2 and receive one unit of asset 1 . Let asset be a riskless 
bond B with annual interest rate r under frequent compounding and face value Bt = 1. We 
assume constant transaction costs Aq for the bond with bid and ask prices as in (j4.5p . Assets 
1 and 2 are two correlated stocks S^ and S"^, denoted in the same currency as the bond, with 
initial stock price for the first stock Sq = 45, volatility o"i = 0.15, constant transaction costs Ai 
and for the second stock Sq = 50, (T2 = 0.2, A2 and correlation p = 0.2 between both stocks. 
The tree is modeled as described in section 14.21 The bid and ask prices are given as in ()4.5p . 
Consider an exchange option with physical delivery. The payoff is given by 

X (Uj) = {Xl ((j) , X2 (Uj) , X3 (W))'^ = (0, /|^a,l>_5,a,2| (W) , -/|^a,l>5,a,2| (w) j . 

The maturity is one year. Table [2] gives the vertices of SHPq{X) in units of (bond, asset 1, 
asset 2)^ and the scalar superhedging prices vtq (X) in units of bond and tt°'{X) in the domestic 
currency for different values for r and A = (Ai, A2, Aa)"^. The recession cone of SHPo{X) is 
equal to Kq generated by the vectors given in (|4.6p . Note that 7r"(X) in the domestic currency 
can be calculated straight forward if Aq = by vr"(X) = itq{X)Bq. If Aq > 0, the scalar 
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superhedging price 7r'^(X) in the domestic currency can be calculated by adding to SHPq{X) 
(as a 3 dimensional object in a four dimensional space, where the cash axis was added) the four 
dimensional cone Kq , which is generated by the bid and ask prices of the bond and the stocks as 
in ()3.9p . and calculating the vertex of the intersection with the cash axis. The reason is that the 
transaction costs for the bond might lead to the effect that there might be cheaper ways to trade 
cash into the set SHPq{X) than to trade the pure cash position {ttq{X))~^Bq — {ttq{X))^ Bq into 
the pure bond position ttqIX), see the fifth example in table [21 where 7r"(X) = 6.988 < 7.011 = 

{Kix)rBS - K(x))-Bg. 

Example 4.4. (outperformance option: superhedging portfolios and strategies) Let us consider 
an outperformance option. Let asset be a riskless cash account with zero interest rates. Assets 
1 and 2 are two correlated stocks S^ and S'^, denoted in the same currency as the cash account, 
with initial stock price for the first stock Sq = 50, volatility ui = 0.15, constant transaction 
costs Ai = 0.2 and for the second stock Sq = 45, a2 = 0.2, A2 = 0.1 and correlation p = 0.2 
between both stocks. The tree is modeled as described in section S21 The bid and ask is given 
as in ([33]). 

We want to calculate the set of superhedging portfolios of an outperformance option with 
physical delivery and strike K = 47. The payoff is given by 

X{u;) = {X^{io),X2{io),Xs{io)f 

-^-^{max{5j'\5j'2}>i^} ('^) ' ^{5j'i>5j'2 and S^-^>K} ('^) ' ^{5j'2>S^'' and Sp^>K} (^) 

The maturity is one year. We will use only a small number of time intervals n = 4 to illustrate 
different possibilities of choosing optimal superhedging strategies as described in section 13.11 
The set of superhedging portfolios SHPo{X) has two vertices 

-27.404 \ / -34.254 

0.514 , 0.567 

0.388 / \ 0.480 

and a recession cone equal to the solvency cone Kq. The scalar superhedging price is 7r"(X) = 
22.624 in the domestic currency and the scalar subhedging price is vr (X) = —8.633. Negative 
bid prices might occur when physical delivery is considered in markets with transaction costs. 
This issue was discussed and resolved in [29], see also remark 3.30 in |31| . 

Let us compute superhedging strategies for a given path using the methods in section 13.11 
We fix a path given by the sequence of indices ji = 1,2,3,3,4 and J2 = 1,1,2,3,4 at times 
t = 0, 1, ...,4 and an initial portfolio vector given by the first vertex xq = (—27.404,0.514,0.388)^. 
Instead of visualizing all possible strategies as a subset of IR , the bi-criteria problems (I3.12P 
and ()3.13p are considered. Plotting the upper images, i.e. the values of the IR^-valued objective 
function and all (component-wise) larger vectors, we can visualize the tradeoffs and select IR^- 
minimal vectors. In particular, all vertices of the upper images are minimal. 

Figure [H illustrates the strategy of following the max-cash approach at all points in time. At 
t = 2, 2.882$ can be withdrawn while still guaranteeing superhedging. However, if one follows 
the min-trading strategy as illustrated in figure [5] and does not take the money out at t = 2, 
one can withdraw even 6.143$ at terminal time i = 4 for this particular path. Thus, by leaving 
money inside, a larger payout at a later date is possible. However, one could also loose the 
money left inside. Figure [6] illustrates this in an example of a mixed strategy: at t = 2 an 
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Table 2: The set of superhedging portfolios of an exchange option with and without transaction 
costs for the bond (see example 4) 
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efficient point in between niax-cash and niin-trading is chosen, and at time t = 3 min-trading is 
chosen which corresponds to a withdrawal of 3.006$, which is 1.536$ less than one could have 
taken out with max cash. However, this money is lost at t = 4. 
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Figure 4: Example 14.41 with max-cash selection. In the first picture we see that no trading (y2- 
axis) is necessary and no withdrawal (yi-axis) is possible at time t = 0. At time t = 1, trading 
is necessary and no withdrawal is possible. We see that for t = 0,1,3,4 exactly one minimal 
point exists, which refers to both a max-cash and a min-trading strategy. At time t = 2 we 
choose point Pi which refers to a max-cash strategy where 2.882$ are withdrawn. A strategy 
that refers to Pi entails a high trading activity. As no further withdrawal is possible, so the sum 
of all withdrawals is 2.882$. 
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Figure 5: Example 14.41 with min-trading selection. The first difference in comparison to figure 
Hlis at time t = 2, where we choose the point P2. The possible withdrawal of 2.882$ is waived 
and a min-trading strategy is chosen. As a result we have different choices at time t = 3 
and a withdrawal of 5.704$ would be possible now. But a better result can be obtained if 
the withdrawal is waived and a min-trading strategy is taken by choosing point P3. At time 
t = T = 4, we have two choices. To obtain the withdrawal only in cash, we choose point P4. The 
sum of all withdrawals is 6.143$. This shows that the min-trading strategy may yield higher 
total withdrawals. 



5 Reformulating and solving the problem 

There are at least four reasons to relate the recursive representation of the set SHPq{X) of 
superhedging portfolios in (|3.2|) and (|3.3|) to linear vector optimization: 

1. The problem can be formulated as a sequence of linear vector optimization problems. 



2. Vectorial duality theory [181I26J in an extended form [16] leads to interesting insights into 
the problem's nature, see also section 17.31 

3. Existing algorithms |31tl32| for the scalar superhedging price as well as their generalization 
to the case of more than two assets can be interpreted in the framework of vectorial duality 
theory, see section [7l 
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Figure 6: Exainple l4.4t loss of money by min-trading strategy. The first difference in comparison 
to figures m and [5] is at time t = 2, where we choose the minimal point P5. The withdrawal at 
i = 2 is 1.222$. At time t = 3 a withdrawal of 3.320$ would be possible connected with trading. 
But we only take the smaller amount 1.7841$ which can be withdrawn without trading (choice 
of point Pq). We see that no further withdrawal is possible at time t = T = 4. This means that 
the total withdrawal is only 3.006$. This means that the money which was waived to withdraw 
at t = 3 by choosing the min-trading strategy has been lost. The max-cash strategy at t = 3 
had result the higher withdrawal 4.542$. 

4. The set of superhedging portfolios can be computed using existing algorithms from linear 
vector optimization (in each iteration step). To this end the algorithms from the literature 
[2l|9l|26] have been extended to more general ordering cones, see [16]. 

We consider a linear optimization problem with a g-dimensional objective function. The 
image space IR'^ is partially ordered by a polyhedral convex cone C C IR'^ that contains no lines 
and has non-empty interior. For y,z € IR"^ we write y <c z, or shortly y<z,ifz — y&C. We 
consider the problem to 



minimize P : IR — t- IR'^ with respect to <c subject to Bx > b, 
where P is linear, i.e. P £ W""^, B G IR'"^'^, h G IR*". The feasible set of ([£]) is denoted 

S:=\xe IR'^I Bx>b\. 

The dual problem to ([P]) is 

maximize D* : IR" x IR'^ — ^ IR'^ with respect to <k over T, 
with (linear) objective function 



(P) 



(D*) 



D* :IR" 



W 



ordering cone 



— ;■ IR"^, D*{u,w) := [wi, ...,Wq^i,b u) , 
is: := IR+ • (0, 0, 



.0,1 



and feasible set 

T := {{u,w) G IR™ X IR5| u > 0, B^u = P'^w, Jw = 1, w e C7+} , 

where c is a fixed vector in int C and C^ := {w G IR"^! \/y £ C : vj^y > O} is the dual cone of 
C. This kind of vectorial duality, called geometric duality, has been introduced in |18j . A more 
detailed discussion is given below in section [6l For the moment we are content with the fact 
that Benson's algorithm [21 [9l |26l [IS] can be used to compute solutions (as defined below) to 
both the primal and the dual problem. The set of superhedging portfolios can be computed by 
the following algorithm. 
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SHP- Algorithm 

Input: 

for all t G {0,1, . . . , T}, for all uj G Qt'- Ktioj) (matrix of generating vectors of Kt{uj)), 
for all UJ S VLt'- X{uj), 
Output: 

for all t G {0, 1, . . . ,T}, for all uj G Jlj: SHPt{X){uj) = {x G IR'^I B^x > 6^}. 



00 


begin 


01 


for all UJ G ^t 


02 


begin 


03 


B^ = Kt{uj)+; 


04 


h^ = (^B^)T . x{uj); 


05 


end; 


06 


for t=T-l downto 


07 


begin 


08 


for all a; G r^t 


09 


begin 


10 


B = {B^\ UJ Gsucc(a;)}; 


11 


6 = {6'^| UJ G succ(tj)}; 


12 


P = LiquidationMap(irt(a;)); 


13 


C = P-Kt{uj); 


14 


compute a solution \^{v},w^) 


15 


B^ = {pTw\...,P^w''f; 


16 


b'^ = {b^u\...,b^u>'f; 


17 


end; 


18 


end; 
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end. 



,(^^^w;'^)} to ([P^: 



The SHP- Algorithm yields the set of superhedging portfolios by computing a sequence of 
linear vector optimization problems together with several elementary operations. To be more 
precise, the algorithm consists of a reformulation of the recursive description of the set of super- 
hedging portfolios given in ()3.2p and ()3.3p . Solving the linear vector optimization problem ()D*P in 
line 14 means to compute an inequality representation of PlSj + C, where S" = {x G IR | Bx > 6} 
with P,B,b,C in each step as in lines 10-13 of the SHP- Algorithm. Under the assumptions of the- 
orem lSTTJ the algorithm computes inequality representations of SHPt{X){uj) for all t G {T, . . . , 0} 
and all uj £ Qf 

Let us discuss the steps of the SHP- Algorithm in detail. In line 03 we need to compute the 
generating vector of the dual cone which is a standard problem of computational geometry. For 
special classes of solvency cones we can exploit an explicit description of the dual cone, see [27] . 
In lines 03 and 04 we store an inequality representation {x G IR | B^x > b'^} of the sets 

SHPt{X){uj) = Xt{w) + Kt{uj) for ah uj G VLt- 

In lines 10 and 11, an inequality representation S = |x G IR | Bx > b\ of 

fl SHPt+i(X){uj), [uj^rtt) 

aJGsucc (oj) 
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is obtained. The set S is understood to be the feasible set of the primal linear vector optimization 
problem ([Pj, where the objective function and the ordering cone are computed in lines 12 and 13. 
The objective function P in line 12 is the liquidation map, which corresponds to an exchange of 
those assets with no transaction costs into a single asset. For instance, if there are no transaction 
costs between assets i and j, the matrix 

P,, = [e\..., e'-i, e^ + 7T'^e^,e'+\ ..., e^-\e^+\ ■■■,e''f (5.1) 

is considered, where vr*-' is the exchange rate, the price of asset j in terms of asset i, as defined 
in section [21 If x is a portfolio with d assets then PjiX is a portfolio with d — 1 assets which is 
obtained from x by exchanging asset j in asset i without transaction costs. The same procedure 
is applied to the new portfolio if there are further pairs of assets having no transaction costs. 
This process determines the liquidation map P for solvency cones as introduced in section [21 
At the end we have transaction costs between any two assets. As a consequence, the cone C 
computed in line 13 (the ordering cone of ([P])) contains no lines. This is important for the usage 
of Benson's algorithm. 

A solution to the vector optimization problem ([P]) can be interpreted as a finite set of efficient 
points (portfolios) and directions that contain complete information to describe the polyhedron 
P[S] + C, which is the set of superhedging portfolios after applying the liquidation map. A 
solution to the dual problem (|D*P provides an inequality representation of the set P[S] + C. An 
inequality representation of SHP{X)t{oj) is then immediately obtained using the liquidation 
map (lines 15 and 16), t G {0, ..., T — 1}. 

Note that we can use Benson's algorithm to compute a solution to ()D*P in line 14. This algo- 
rithm yields simultaneously a solution to ([P]) . This means that we also obtain a representation of 
the superhedging portfolios SHPt{X){Lo) by finitely many efficient portfolios |x^, . . . , x*| C IR 
and efficient directions |x^, . . . x*} C ]R \ {0}, that is, a representation of the form 

SHPt{X){Lo) =co {xS...,x'}+cone {x\ . . .x^} + Kt{Lo). 

A definition of what is meant by efficient is given in the next section. 

It should be noted that line 14 could be replaced by methods from computational geometry. 
Essentially one can use vertex enumeration methods for (unbounded) polyhedra. Benson's 
algorithm also involves a vertex enumeration method and additionally several LPs to be solved. 
Beyond the fact that Benson's algorithm can be easily applied, we obtain interesting insights 
into the nature of set-valued superhedging portfolios by relating them to vector optimization. 
Therefore, we continue this discussion after introducing same basic notions and results from 
linear vector optimization. 

Note that the output of the SHP- Algorithm can be used to compute superhedging strategies 
as discussed in section [3Tl where we assumed that inequality representations of SHPt{X) are 
known. 

6 Some basics on linear vector optimization 

It should be noted first that there are different approaches to solution concepts and duality for 
linear vector optimization problems in the literature, see e.g. [201 [3 SI [26] and the references 
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therein. Here we only present geometric duality [181 [Ml US] and solution concepts as introduced 
in [26]. 

Let us first recall some facts on polyhedra. A polyhedron (or polyhedral set) A in IR'^ is 
defined to be the intersection of finitely many halfspaces, that is 



A=f]{yeW\iz'fy>j,} (zVO). 



i=l 



Every non-empty polyhedron A in IR'^ can be expressed by finitely many points x^, . . . ,x^ € IR^ 
and directions A:\ . . . , A:* S IR^ \ {0}: 



t 



A= <]^AiX^ + ^;UjA:^' A,>0(i = l,...,s), X]^* = l' /^j > (j = 1, . . . , t) 



=1 



> . 



This can be also written as 

^ = CO {x\... jx"} +cone {k^,... ,k*} . 

Note that we set cone0 = {0}, which applies in the bounded case where no directions occur. The 
set Aoo '■= cone jA;^, . . . , A;*} is the recession cone of A. As usual in computational geometry, a 
finite set of halfspaces defining a polyhedron A is called //-representation (or inequality repre- 
sentation) of A, whereas a finite set of points and directions defining A is called ^/-representation 
of ^. A bounded polyhedron is called poly tope. 

Let us consider problem ([P]) as introduced in section \E\ with objective matrix P G IR'^^ , 
polyhedral ordering cone C C IR'^ that contains no lines and has non-empty interior, and feasible 
set 5 = {x E IR*^! Bx > b}. A solution to ([P|) (as defined in [26j) is a finite set of points and 
directions in S being efficient for ([P]) (as defined below). Moreover, a solution to (JP]) is required 
to provide a ^/-representation of the polyhedron V = P[S] + C . 

To be more precise, a point x E S* is said to be efficient for ([P]) if there is no z E S" with 
Pz <c Px and Pz ^ Px. The set of efficient points of (JPJ) is denoted by Eff([P|). Efficient 
directions can be introduced by the homogeneous problem, 

minimize P : IR — )■ IR"^ with respect to <c subject to Bx > 0. (P'*) 



Nonzero efficient points of (jP |) are called efficient directions of (jP|) 



A solution t o (|Pl) is a non-empty finite subset S C EffdP]) together with a (possibly empty) 



subset S^ C EffdO such that 

CO P[S] + cone P[S^] + C = P[S] + C. 

The latter condition means that a ^/-representation of the set V := P[S] + C is given by a 
solution. The set V is called upper image of (jP]) . 

Likewise a feasible point {u, w) £ T is said to be efficient for (|D*j) (as introduced in section 
[5|) if there is no (n, w) gT such that D*{u, w) >k D*{u, w). Efficient directions do not occur in 
the dual problem. The set of efficient points for (ID*j) is denoted by Eff (lD*j) . A solution to (ID*|) 
is a non-empty finite subset T C Eff (ID*j) such that 



CO D* [T] - K = co D* [T] - K. 
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This means that a solution can be used to get a ^/-representation of the set D* := D*[T] — K, 
which is cahed the lower image of ()D*p . 

To give a first flavor of the duality relation between ([P]) and ()D*p . let us mention here that 
a solution to (JP]) also provides an i^-representation of V* = D* [T] — K and a solution to (|D*|) 
provides an //-representation oi V = P[S] + C , which is exactly what we use in the above 
algorithm. 

The duality relation between V and V* is similar to the classical idea of duality of polytopes. 
Two polytopes A and A* in IR*^ are said to be dual to each other if there exists a one-to-one 
mapping ^ between the set of all faces of A and the set of all faces of A* such that ^ is inclusion 
reversing, i.e. faces F^ and F^ of A satisfy F^ C F^ if and only if the faces ^(F^) and ^(F^) 
of A* satisfy ^(F^) D ^(F^), see e.g. [12]. 

A duality mapping ^ for the vector optimization problems ([P]) and (jD*|) is now introduced 
in several steps. We assume throughout that there exists a vector 

c G int C such that Cq = 1. (6-1) 

For the parameter c, we consider the following bi-affine function 

q-l / q-l \ 

y. : M^ X M^ ^ IR, if{y,y*) ■.= Y,yiyl ^ y^ 1 - ^ Qy* - y\. 

i=\ \ i=\ I 

This coupling function has been introduced in [18] for the choice c = (1, ..., 1) . 

The coupling function ip is used to define the following two injective hyperplane-valued maps 

R:^^^ ^\ H{y*) := {y G IR"! ^(y, y*) = 0} , 

H*:W^ IR^ H*{y) := {y* G IR^| ^{y,y*) = 0} . 
The mapping H is used to define the duality function 

^.2iR'^2'^', ^(F*) := Pi H{y*)nV, 

yaF* 

where 2^^ denotes the power set of IR"^. 
Weak duality reads as follows. 

Theorem 6.1 (|181l26j). The following implication is true: 

[yeV^y*eV*] =^ (^(y,y*)>0. 

A motivation can be given by the following geometric interpretation. Weak duality implies 
the inclusions 

V* C {y* G IR^I Vy G P : if{y, y*) > 0} and P C {y g IR«| Vy* G V* : ip{y, y*) > 0} , 

whereas the following strong duality theorem yields even equality. 

Theorem 6.2 ([TBI [26]). Let the feasible sets S of ^ and T of ([E]) be non-empty. Then 

[Vy*GP*:v?(y,y*)>0] ^ yeV, 
[VyGP:<^(y,y*)>0] =^ y* e V* . 
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Geometric duality is a further duality statement that takes into account the facial structure 
of the sets P and P*. 



Theorem 6.3 ( [1S\ I17j). ^ is an inclusion reversing one-to-one map between the set of all 
K -maximal proper faces ofT>* and the set of all proper faces ofV. The inverse map is given by 

^-\F)= f]H*{y)nV*. 

y&F 

Moreover, if F* is a K -maximal proper face ofD*, then 

dimF* + dim^(F*) = q - I. 

As shown in [26], the non-iT-niaxinial proper facets (i.e. faces of dimension q — 1) of D* 
correspond to the extreme directions of P by a similar relation, where the coupling function has 
to be replaced by ip : IR'' x IR'' — )• IR, ip{y, y*) := ip{y, y*) + y*. Let us illustrate the main ideas 
by two examples. 
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Figure 7: Illustration of example 16.41 The left two pictures show the feasible set 5, its image 
P\S\ and the ordering cone C . The right two pictures show the upper image V = P[S] + C of 
the primal problem ([P]) and the lower image P* = Z)*[r] — i^ of the dual problem (|D*[) and 
display the idea of geometric duality: The three facets (edges in the current 2-dimensional case) 
of V correspond (via ^) to the three vertices of D* . We distinguish ET-maximal (plain lines) 
and non- if -maximal facets (dashed lines) of T>*. The two if -maximal facets correspond to the 
two vertices of V and the two non- if -maximal facets correspond to the two extreme directions 
of V (via ^ defined likewise to ^ using ip). 



Example 6.4. Consider problem (JP]) with the data 



P 



B 



and the ordering cone 



C 



cone 



/ 2 


1\ 


1 


2 


1 





\o 


1/ 


-3\ 


( ^ 


1 / 


' 2 



6 




The pair {S,S^) with S = {(2,2)^, (6,0)"^} and S^ = {(l,0)'^}_is a solution to (|P]). It contains 
the information to construct the upper image V via V = co P[S] -|- coneP[S''^] -|- C, see figure [7] 
for an illustration. 
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Choosing c = (0,1) G int C, we obtain the geometric dual problem (|D*P with the objective 
function 

D*{ui,U2,U3,U4,Wi,W2) = {wi,6ui +6^2)^ 

and the feasible set T C IR x IR which is given by the equalities 

2ui + U2 + U3 — Wl — -0)2 = 

Ui + 2U2 + U4 + Wi — -0)2 = 

W2 = I 

and the inequalities 

ui > 

U2 > —3wi + ^2 > 

1*3 > lUi + 2^2 > 

ti4 > 

The set 

t= i(0,0,0,2,-l,l)'',(o,|o,0,-i,l') ,[1,0,0,0,^,1 



provides a solution to (|D*p . It contains the required information to construct the lower image 
V* := D*[T] - K via V* = D*[f] - K. The vertices of V* are 

O-Ifl = {(-l,0f,(-i.4)',(i,;' 

and together with the extreme direction (0, 1) of the cone K we have a IZ-representation of V* . 
Using the coupling function (p we immediately obtain an //-representation of the upper image 
V. For example, taking y*^ = (—1,0)"^ G D*\T], the inequality ip{y,y*^) > yields —2/1 + 2/2 > 0. 
An //-representation of V = P[S] + C is obtained in this way as 

V = LeJR^ - 2/1 + y2 > 0, --2/1 + 2/2 > 4, -2/1 + 2/2 > 4 1 . 

Likewise we obtain an //-representation of V* from a solution to (P). For x^ := (2,2)-^ G S 
and y^ := Px^ = (0,4)-^ the inequality ip{y^,y*) > yields 2/2 < 4 and for x^ := (6,0)^ G S 
and 2/^ := Px"^ = (6,6)"^ we get dy^ — 2/2 ^ ~6. Further inequalities are obtained from S^ 
and from the extreme directions of the ordering cone C. To this end, <p has to be replaced 
by if, which means that the 2/2-component has to be deleted. Taking x^ = (1,0)"^ G S^ and 
yi := Px^ = (1, 1)^, we get yl > -1. Taking y"^ = (-3, 1)^ and f = (1, 2)"^ we get 2/t < i and 
2/1 > —2, where the latter inequality is redundant. Thus 



P* = <! 2/* G ]R^ 



2/2 < 4, 62/1 - 2/2 > -6, -1 < 2/1 < 3 



see figure [71 
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Figure 8: The upper and lower images of Example 16. 51 The dark gray facets, the plain edges 
and all vertices of P* are i^-maximal and the light gray facets and the dashed edges of T>* are 
not. 



Example 6.5. Consider problem (JP]) with the data 
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and the dual problem (ID*P for the parameter c = (1, 1, 1)^. The upper image V has 6 vertices, 18 
edges and 13 facets. The lower image D* has 13 vertices, 18 iT-maximal edges and 6 X-maximal 
facets, see figure [8] for an illustration. 

Note that in [181 126] only the case c = (1, ..., 1)^ and C = IR!?^ is considered. The proofs of 
the general case are straightforward and will be published in a forthcoming paper [16]. Moreover, 
Theorem l6.3l (for c as in ()6.ip ) is a special case of geometric duality for convex vector optimization 
problems, see example 3 in |17j . 



7 Scalar superhedging price 

The set of superhedging portfolios plays an important role when one is actually interested in 
carrying out a strategy starting from an initial portfolio vector that can contain more assets 
than just cash. Clearly, in this case, it is not enough to know the scalar superhedging price. 

On the other hand, if one is interested solely in price bounds for a claim in a certain currency, 
the scalar superhedging and subhedging prices 7r'^(X) and '7r^(X) give exactly this information. 
Prom robust no arbitrage it follows that the market bid and ask prices p {X) and p°'{X) of a 
claim X have to satisfy 



p\X) < TT^iX), p^iX) > tt\X) and p\X) < p^^X), 



(7.1) 
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where p {X),p"'{X),'K {X),'k'^{X) are all denoted in the same currency. Furthermore, the in- 
equalities 

p\X) > vr^(X), and p"(X) < ^"(X) 

are reasonable: Obviously, one would rather superreplicate X with tt"'{X) than to buy it at a 
higher price if p'^{X) > 7r"(X), but this last inequality would not create robust arbitrage as long 
as the inequalities in ()7.ip are satisfied. In total, one obtains price bounds 

tt\X) < p\X) < p^iX) < TT^iX). 

In the following, we will discuss how vr''(X) and 7r"(X) can be calculated and how the obtained 
algorithm is related to the previously studied SHP- Algorithm. In section 17. H theorem I7.H a 
dual representation of the scalar superhedging price is given. The result is the d-dimensional 
version of Jouini, Kallal |21] and can be obtained by scalarizing the dual representation of the set 
of superhedging portfolios. The dual representation of the scalar superhedging price allows to 
deduce dynamic programming equations (corollary 17. 2p that allow to implement and efficiently 
calculate the scalar superhedging price of a claim (algorithm in corollary 17.31 in section 17. 2p . 
In section 17.31 the relation between the scalar algorithm and the SHP- Algorithm via geometric 
duality is discussed. One obtains that the calculation of vr"(X) reveals also the set SHPq{X) 
if a certain mapping is applied to a function appearing in the penultimate step of the scalar 
algorithm (lemma l7.5p . 

7.1 Scalarization of the dual representation 

In the previous sections, the set SHPq{X) of all initial portfolio vectors that allow to superhedge 
a claim X G L^{Tt,^ ) was studied. If one is interested in the calculation of price bounds, it 
is helpful to study the smallest superhedging prices (and the largest subhedging prices) in the 
currencies or numeraire of interest. To do so, consider the one dimensional subspace M (see also 
(j4.ip ) given by M = {se*,s E IR}, where asset number i is the chosen numeraire. One might 
repeat this procedure for different currencies/numeraires if one is interested in price bounds of 
the portfolio X in different currencies/numeraires. We focus on those superhedging elements 
that lie in M, i.e. we consider SHPq{X). Prom ()2.3p and theorem 12. H we obtain 

SHP^{X) = {xo£M: X exo + Ar} (7.2) 

Pi {E^ [X] + G (w)) n M (7.3) 

{iQ,w)eW^: EQ[X-]<oo} 

= {xo e M: VZ G Z with E[{X^ Zt)'] < oo : E[X^ Zt] < x^ Zq} . (7.4) 

To obtain the smallest superhedging price, we apply the scalarization procedure introduced in 
[13] to the function R{X) = SHPq^{X). That is, we consider the extended real- valued function 
m,v ■■ L^di^T, IR'^) ^ IR U {±00} given by 

<fR,v (X) = inf v'^u 
u&R{X) 



for V G {Kq )^ (see also section 5.1 in [15j). {Kq )'^ denotes the positive dual cone of the cone 

{K^)+ = {v £ M:yu£ K^^: v^u > O} C M. 



K^^ = KonM inM. Thus 
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We will apply the scalarization to the dual representation of SHPq^ given in (j7.3p , respectively 
(j7.4p . Dynamic programming equations can be obtained. For finite probability spaces this leads 
to an algorithm that goes backwards in the event tree. Recall that in our case M = {se*, s S IR}, 
where asset number i G {1, ..., d} is the asset of interest, e.g. the USD cash account, or a bond. 
Thus, we are scalarizing with respect to v = e* G {Kq^)~^ = {se*,s G IR+}. The calculation of 
the smallest superhedging price in the asset of interest, that is the calculation of vrf (X) leads 
to the a generalization of the well known Jouini, Kallal [2T] representation to the d asset case. 
For simplicity we assume that the solvency cone Kt contains no lines. Then, the solvency cone 
Kt{(jj) is spanned by the vectors vr^-'e* — e^ , 1 < i,j < d, see section [2j The general case can be 
derived using the liquidation map in (j5.ip . The following theorem gives a dual representation 
of the scalar superhedging price in the spirit of Jouini, Kallal |21] in the d asset case. 

Theorem 7.1. Under the robust no arbitrage condition (NA'^), the scalar superhedging price 
7r?(X) in units of asset i G {1, ...,d} is given by 

<(X) = sup E^[X^St], (7.5) 

where Q* is the set of all processes {St)J=Q with SI = 1, S^ < ii^^Sl for all 1 < j,k < d and 
their equivalent martingale measures Q for which E'^KX St)~] < oo. 

Proof. From the scalarization of (j7.4p with respect to w = e* G (i^g'^)+ one obtains 
<(X) = inf v'^u 

u<^SHP^''{X) 
= min{tG]R: sup E[X'^^]<t}. (7.6) 

{Z &Z:E[{XT Zt)- \<oo} ^0 

Note that for every Z G Z, one can define the corresponding frictionless price of the d assets 
expressed in asset i as 

yl 72 yd 

5^ = (fl,fL,...,fL), (7.7) 

i.e., SI = l,and obtains an equivalent martingale measure Q of the process (5j)^q via 

dQ _ Z'j. 
dP~~Z^' 

Note that Zl > 0, t G {0, ...,T} is ensured by assumption 12.11 The set Z of all consistent price 
systems Z is one-to-one (up to a multiplicative factor for Z) to the set of processes (6'i)^Q with 
SI = 1, S^ < TT^ Si for all 1 < j,k < d and their equivalent martingale measures Q. This follows 
from the observation that St is defined by Z G Z with Z G K^\ {0} via (j7.7p and Kt is spanned 
by the vectors vr'-'e* — e^ , 1 < i,j < d. Furthermore, Z G Z is a martingale under P, which 
corresponds to Q being a martingale measure for (5i)^g, see [SQ p. 24/25 for details. Since for 
every Z G Z with E[{X'^ Zt)'] < oo 

E[X^^]=E[X^^^]=E<^[X^Sr], 



we can rewrite (I7.6P as in (|7.5p . D 
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The supremum in (|7.5p is attained if we replace Q' by the enlarged set Q , that contains 
martingale measures that are not necessarily equivalent to P . We will write dynamic program- 
ming equations for problem (j7.5p that will allow to implement and efficiently calculate the scalar 
superhedging price of a claim X G L^(J^r,]R ). 

To do so, let us define the the sets of one-step transition densities as in [7] 

Vt := le G Li(j;,]R+) : Et-i^ = 1}, t = 1, ...,T. 

Let Co = 1- Every sequence ((S't,Ct))?lo ^ith Sf = 1, S^ < tt^^sI for all I < j,k < d for 
t = 0, ...,T and £,t £ T^t, Et-i[^tSt\ = St-i for t = l,...,r defines an element {St,Q) in Q by 
setting 

On the other hand, every element {St, Q) ^ Q induces a sequence with the above properties by 
setting 



I 1 ^T1 J F . , 

Up 
Let us denote for t = 1, ..., T 



1^ on{i?i_i[^]>0} 

dP' 

1 on{^,_i[f]=0}. 



^j(5i_i) = {{St,Ct) : S,* = l,5f < 7T^'s{,l<j,k< d, 6 G A, i?t-i[65t] = 5t-i} • 

From the considerations above, we obtain the following. 

Corollary 7.2. Assume the robust no arbitrage condition (NA^). The scalar superhedging price 
vrf (X) given in (I7.5P can be written as a sequence of nested optimization problems. Let the value 
function at time T — 1 be 

Vt-i{St-i) = max Et^i[CtX^ St]- (7.8) 

For t G {T — 2, ..., 0} we define the value function 

Vt{St)= max Et[^t+^Vt+i{St+i)]. (7.9) 

Then 

<(X) = max Vo{So). (7.10) 

So&'SR'^:Si^=l,S^<TTJ'=S^Q,l<j,k<d 

Proof. This follows from the one-to-one correspondence between the set Q and {{{St,Ct))t=o '■ 
S'o = l,S^ < TT^^Si, l<j,k<d, (Si, 6) G Ai{St-i),t = 1, ..., T} and the tower property. D 

In the following section an algorithm is described that is based on this recursive equations. 
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7.2 Algorithm for the scalar superhedging price 

Let us assume a finite filtered probabifity space {^,T, {Tt)f=o,P) as in section [31 We will split 
each of the iteration steps in (j7.9p into two substeps. In the first steps, one incorporates only the 
constraint Sl_^_l = 1, S^^^ < tt^^S-J-^i for 1 < j,k < d, the second steps coincides with actually 
solving (|7.9p . For d = 2, this algorithm coincides with algorithm 4.1 in |32j . and algorithm 3.15 
in ED. 



Corollary 7.3. (Algorithm scalar superhedging price under assumption (NA)) 
1. 

{ X^St : Sip = l,Sl^< TT^'^S^r, l<j,k<d 
Vt{St) = Vt{St) = < (7.11) 

I — oo : else. 

This defines a function V^ : IR — )• IRU {— oo} at each node uj G ilr- 
2. For t G {T — 1, ..., 0} and nodes u £ i^t 

V^^iSt) = cap{V,%^{St+i) : u G succ (u)}. (7.12) 



Vt{S, 



VtiSt) : Sl = l,S'^<7T^''Sll<j,k<d 
— cx) : else. 



(7.13) 



3. The scalar superhedging price of X £ L^(J^t',IR ) is given by 

<(X)= max Vb(5o). 

The concave cap function cap {/i, ..., fm} : IR — )■ IR U {— oo} of a finite number of concave 
functions /i : IR — )■ IRU {— oo}, i = 1, ...,m is defined by its hypograph via 

hypo(cap{/i,...,/m}) =co I (J hypo /i I , (7.14) 

compare [31], p. 135. The closure can be omitted if /,, i = l,...,m are polyhedral, i.e., the 
hypographs of /j are polyhedral convex sets. 

Proof, (proof of corollary 17. 3p It follows from the representation of the cap function as an 
optimization problem as in lemma A. 4 in [3T] that the two substeps (|7.13p (respectively (|7.1ip 
for t = T) and (j7.12p coincide with optimization problem (j7.9p . D 

Remark 7.4. Obviously, the scalar superhedging price at time t and node w G Oj is given by 

«)t(X)(a;) = max y-(SO = max , V^St). 

Stew^ St£W^:Si=i,si^<TTl'^{uj)sl,i<j,k<d 

In section [7^3] below, we will discuss how the set of superhedging portfolios SHPt{X) at time t 
and node a; G ^2^ can be obtained from V^{St) via geometric duality. 

Equation ()7.14p already indicates a relationship between the cap function and set-operations. 
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7.3 Interpretation in terms of vector optimization and recovering the set of 
superhedging portfolios 

We show in this section that the algorithm for the scalar superhedging price is closely related 
to the SHP- Algorithm which computes the set of superhedging portfolios. This means that 
SHPq{X) is also obtained by the scalar algorithm and, depending on how the scalar algorithm 
is realized, it practically coincides with a special case of the SHP- Algorithm. This relationship 
is established using the linear vector optimization reformulation of the recursive representation 
in theorem 13.11 and geometric duality, see sections [5] and [6l 

Lemma 7.5. Let the assumptions of theorem \3.1\ be satisfied and let V^^^St) be as in corol- 
laryH^ For t = 0, ...,T and uj G ^t, 



SHPt{X){u) = {x e IR'^I {Stfx > Vt^iSt), (SuVt'^iSt)) vertices o/hypoT4}, (7.15) 

Proof. We start with a reformulation of the scalar algorithm. As it can be seen from corollary 17.21 
and 17.31 using the first three substeps, one obtains for uj £ ^t-i 



VjLi(5't-i) = < 



rnax ^ S.'^X{ofS^ if S'^_^ = 1, S^_-^ < 4-i(^)'5^-i, l<j,k<d 

iIiGsucc (a;) 

— oo otherwise 

(7.16) 

where the maximum is taken subject to the constraints 

{s^y = i,{sn''<7T?r\o){sn\i<j,k<d, r>o, E ^" = 1' E r5^ = ^T-i 

[DSsucc (a;) dJSsucc (lj) 

For f = T — 2, . . . , and cj G $7^ one has 

Vr{St) = l («"'^")^gsuccH (717) 

— oo otherwise, 

where the constraints for the maximum are 

r>o, ^ r = i, E ^"3^ = St. 

oJSsucc (oj) iZ'Ssucc {uj) 



Finally, we have 



<(X) = max Vo{So). (718) 

5oGlR'* 



The constraints SI = 1, S^ < irl {(^)Sl, I < j,k < d can be equivalently expressed as S'j = 1, 
St S K^. The constraints in (|7.16p concerning (^'^, 5"^) can be replaced by ^'^5'^ = K^{uj)u^ , 

u^ > Oi T.6J&ucc{u>)^T(^h^ = St-1 using the fact that we set Sfp_^ = 1. From (fflBD we 
obtain for to G ^t-i, 

max Y^ X{uj)^k^{u))u^ ii St-i e K^_^{uj), 5^_i = 1 

V^_i{St-i) = \ "" Q£svicc{ui) (7-19) 

-oo otherwise, 
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where the maximum is taken subject to 

(Zigsucc (uj) 

Setting B^ = {k+{oj)Y, b"^ = {K+{Cj)Y X{u}), w = (S^_i, ..., S^"_\, 5^+_\, ..., 5^_^, 5^_i) and 
omitting the i-th variable S^_^ of Vt-i (which does not changes the problem), we see that 
hypoVr-i is nothing else than the lower image V* of a dual vector optimization problem (|D*P 
where the parameter vector c, which has to be an interior point of the ordering cone, is chosen 
to be the d-th unit vector, i.e. c = e . Note that the i-th component of St-i, corresponds to 
the d-th component of w. This means that c = e'^ corresponds to the numeraire-component. 
The corresponding primal problem ([P]) is to 

minimize id-.M"^ ^ JR'^ w.r.t. <KT-iiuj) s-t. B^x >b^,cd£ succ (a;). (7.20) 

We know by theorem [XT] that SHPt{X){0) = {x £ M'^\ B'^x > h"^} and the upper image 
of problem ()7.20p is "P = SHPt-i{X){uj). Geometric duality (theorem 16. 3p yields that an 
iJ-representation of V is given by the vertices of T>* , that is, 

SHPT-iiX)ioj) = {xe m!'\ {Sr-ifx > V^_^{St-i), (7.21) 

(5t-i, V^7v,i('S't-i)) vertices of hypoVr_i| 

which can be expressed by a matrix B^ and a vector b^, i.e. 

SHPt-i{X){uj) = |x E IR'^I B'^x > 6"^} . (7.22) 

At time t G {T — 2, . . . , 0} we use elements (St+i, Vt+i{St+i)) of the graph of 14+i, which can 
be expressed as a convex combination of vertices of hypoVj+i. The coefficients of the convex 
combinations can be interpreted as variables u of the (geometric) dual ()D*P of the linear vector 
optimization problem 

minimize id ■.M'^ -^ WC^ w.r.t. <Kt(oj} s.t. B'^x >b'^,CJ£ succ (cj). 

Thus (j7.17p can be reformulated as 

f max V {b^)'^u^ if5jGi^+(w), 5| = 1 

yt^{St)={ "" <Desucc(a;) (7.23) 

(^ — oo otherwise, 

where the maximum is taken subject to the constraints 

ti"^ > 0, Y, {B^fu^^ = St- 

aJGsucc (a;) 

We see that hypo V^ is again the lower image D* of the dual vector optimization problem ()D*P 
for c = e"^. Likewise to above, using theorem 13. II and geometric duality (theorem 16. 3p . we obtain 
()7.15p . which completes the proof. D 
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Remark 7.6. Note that (|7.16p and (jT.lTp are parametric linear optimization problems. For 
every choice of the parameter St, a linear program has to be solved. It is well known that 
parametric linear problems of the present type correspond to linear vector optimization problems, 
see e.g. [TT] . 

The above considerations show that the scalar algorithm is closely related to the SHP- 
Algorithm if the last step (jT.lSp is omitted. The SHP- Algorithm computes SHPq{X) by a very 
similar construction if the parameter c = e is chosen and the input data are transformed such 
that the last component refers to the numeraire asset. In contrast to the SHP-Algorithm, the 
liquidation map does not occur in lemma ED and its proof. The reason is that geometric duality 
is not restricted to polyhedral sets that contain no lines, but Benson's algorithm is. On the 
other hand it remains open in lemma [731 how the vertices of hypo 14 are computed. 

Remark 7.7. A relation between the function whose epigraph is the set SHPt{X) and the 
function Vt could also be deduced via conjugation (Legendre-Fenchel transform) in the spirit of 
section 4.2 in |33j . A similar construction has been used to prove a geometric duality theorem 
for convex vector optimization problems [17j . 

8 Conclusion 

In this paper, an algorithm in the numeraire-free approach of Kabanov j22j . Schachermayer [35] 
is presented to calculate the set of all superhedging portfolios in a d asset market with pro- 
portional transaction costs. Furthermore, an algorithm to calculate the superhedging strategy 
when starting from an initial portfolio vector in the set of superhedging portfolios is given. The 
d asset case also allows to consider basket options. Examples are given. 

It is shown that the superhedging problem in markets with transaction costs leads to a 
sequence of linear vector optimization problems. This generalizes the well known fact that in 
frictionless markets, superhedging leads to a sequence of linear optimization problems. 

The above algorithms are the first of its kind. Previous papers, mostly in d = 2, focused 
on the representation and calculation of the scalar superhedging price, which is important to 
calculate price bounds, but does not give sufficient information if one wants to carry out a super- 
hedging strategy starting from portfolio vectors. We extend those results to the d dimensional 
case and show how scalar algorithm like [32^ [3T] or [3] and representation results like |21j are 
related to the set-valued results presented here. It turns out that the hypograph of a function 
appearing in in the penultimate step of the scalar algorithm in [32^ l3T] is the geometric dual (in 
the sense of vector optimization duality) of the set SHPq{X). This allows to recover the set 
SHPq{X), not from the scalar superhedging price 7r"(X) itself, but from the penultimate result 
in the scalar algorithm. 

In [33] and the recent work [34] set- valued constructions to calculate the set SHPq{X) are as 
well presented, but without the connection to vector optimization and corresponding algorithms, 
which leave the problem open how to implement these set-valued constructions for d > 2. The 
approach taken here, answers this question. 
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